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We investigate the phenomenon of particle creation of the massless scalar field in the model of 
spacetime in which, depending on the model's parameter a , the chronology horizon could be 
formed. The model represents a two-dimensional curved spacetime with the topology R 1 x S 1 
which is asymptotically fiat in the past and in the future. The spacetime is globally hyperbolic and 
has no causal pathologies if a a < 1, and closed timelike curves appear in the spacetime if a > 1. 
We obtain the spectrum of created particles in the case a < 1. In the limit a — > 1 this spectrum 
gives the number of particles created into mode n near the chronology horizon. The main result 
we have obtained is that the number of scalar particles created into each mode as well as the full 
number of particles remain finite at the moment of forming of the chronology horizon. 

PACS numbers: 0420G, 0462 
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I. INTRODUCTION 

In 1988 Morris, Thorne, and Yurtsever had demon- 
strated how one could manufacture closed timelike 
curves in a spacetime containing two relatively moving 
traversable wormholes. That work kindled considerable 
interest in the question whether it is possible, in princi- 
ple, to construct a "time machine" — i.e. whether, by 
performing operations in a bounded region of an initially 
"ordinary" spacetime, it is possible to bring about a "fu- 
ture" in which there will be closed timclikc curves. So 
far, in spite of many- year attempts to answer this ques- 
tion, there does not exist full and clear understanding 
of the problem. On the first stage of investigations, it 
seemed that the quantum field theory could give a mech- 
anism which would be able to protect chronology. Hawk- 
ing has suggested the chronology protection conjecture 
which states that the laws of physics will always prevent 
the formation of closed timelike curves |^,^|. His argu- 
ments were based on the supposition that the renormal- 
ized vacuum expectation values of a stress-energy tensor 
of a quantum field must diverge at the chronology hori- 
zon which separates a region with closed timelike curves 
from a region without them. Various attempts at proving 
Hawking's conjecture have been made culminating in 
singularity theorems of Kay, Radzikowski, and Wald jfj). 

However, in recent time, a number of examples of con- 
figurations with the bounded renormalized stress-energy 
tensor near the chronology horizon has been given M Ell . 
For instance, in our works (^|J9) it was shown that, in the 
case of automorphic fields, there exists a specific choice of 
quantum state for which the renormalized stress-energy 
tensor vanishes at the chronology horizon (in fact it is 
zero in the whole spacetime); Krasnikov (?]] exhibited sev- 
eral 2D models, and Visser pi presented a 4D model of 



spacetime - the "Roman ring" of traversable wormholes 
- for which the vacuum polarization can be made arbi- 
trarily small all the way to the chronology horizon; in ad- 
dition, Li |li"| ] showed that the renormalized stress-energy 
tensor may be smoothed out by introducing absorption 
material, such that the spacetime with a time machine 
may be stable against vacuum fluctuations. 

These examples indicate that it is impossible to prove 
the chronology protection conjecture taking into account 
only effects of the vacuum polarization. Recently Visser 
]l2| , ^0]| has given arguments that the solving of the prob- 
lem of chronology protection is impossible within the con- 
text of semi-classical theory of gravity and requires a fully 
developed theory of quantum gravity. 

Nevertheless in this work we remain in the framework 
of semi-classical quantum gravity. We shall discuss such 
the aspect of the quantum field theory in spacetimes with 
the time machine which was so far not investigated. Up 
to this time, various aspects of vacuum polarization near 
the chronology horizon have mainly been investigating. 
Here we consider dynamics of the chronology horizon 
forming. As is known, time-dependent processes in the 
presence of quantum fields are accompanied by a pair 
creation. So one may expect that particles of quantized 
fields will be created in the process of the chronology 
horizon forming. The aim of this work is to determine 
the spectrum of particles and try to answer the question: 
Can the particle creation stop the chronology horizon 
forming? 

We use the units c — G = Ti = 1 throughout the paper. 



II. A MODEL OF SPACETIME 

Here we shall discuss a model of spacetime in which 
the chronology horizon is being formed. Let us consider 
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a strip {77 G (—00, +00), £ s [0, L]} on the ry-£ plane and 
assume that the points on the bounds 7~: £ = and 
7 + : £ = L are to be identified: (rj, 0) = (77, L). After this 
procedure we obtain a manifold M with a topology of 
cylinder: R 1 x S 1 . Introduce on this manifold the metric 



ds 2 



dn 2 + 2a{rj)dr}d£, - (1 - a 2 {r/))d£ 2 



(1) 



where a(r)) is a monotonically increasing function of 77 
which has the following asymptotic behavior: 



a(rj) 







if -q 
if 77 



-00, 
-00, 



(2) 



where a is some constant. Further we shall call the man- 
ifold M as the spacetime M with the metric ([!]). Note 
that the metric 's coefficients in Eq.(0) do not depend on 
£, so they themself and their derivatives are taking the 
same values at the points (r/,0) and (77, L). Hence, the 
internal metrics and external curvatures at both lines 7" 
and 7 + are identical. This guarantees the regularity of 
the spacetime M. 

The metric ([!]) describes a curved spacetime which 
is asymptotically flat in the past, 77 — > — 00 (the "in- 
region"), and in the future, 77 — > 00 (the "out-region''). 
Really, as follows from if 77 — ► —00, the metric (|l|) 
takes exactly the Minkowski form 



ds 2 



drf - d£; 



In the future, n — > 00, the metric (|l|) reads 

ds 2 = dn 2 + 2a dr)d£ - (1 - a 2 )d£ 2 , 
and takes the Minkowski form 

ds 2 = dt 2 - dx 2 
in new 'Minkowski' coordinates 

t = rj + a Q £, x = £. 



(3) 



(4) 



(5) 



(6) 



Note that the spacetime M could be considered as the 
factor space: M = M/1Z. Here M is an universal cov- 
ering spacetime for M, in our case it is the whole plane 
(77, £), and 7^ is the equivalence relation: 



(7) 



Consider now the causal structure of the spacetime M. 
With this aim we have to define null (lightlike) curves 
which form a light cone at a point. The equations for 
null curves could be found from the condition ds 2 = 0. 
There exist exactly two null curves which go via each 
point of the spacetime and their equations read 



dfj 



1 + a(fj) 



const, 



£ — I = const. 

1 - a(rfj 



(8) 



As follows from the asymptotical properties (g) of a(rj) 
in the in-region the equations (@) take the form: 



77 + £ = const, rj — £ = const, 
and in the out-region they are 



(9) 



7/ + (1 + a )£ = const, rj — (1 — a )£ = const. (10) 

Analysing of the equations (||) together with their asymp- 
totical forms ([)]) and (10|) reveals the following quali- 
tative picture: The light future cone in the past (see 
Eqs.(^)) has the angle a\ — 90° and has no inclination 
(i.e., the angle between the cone's axis and the 77-axis is 
zero); then, the cone is enlarging and inclining so that 
in the future (see Eqs.(|l0|)) its angle becomes equal to 
Q'i = 7r — arctan(l + a a ) — arctan(l — a ) and the inclina- 
tion's angle becomes equal to a-i — i[arctan(l + a a ) — 
arctan (1 — a D )]. There are two qualitatively different 
cases: (i) a Q < 1 and (ii) a a > 1. In the first case, 
a a < 1, this rotation of the cone on the 77-^ plane does 
not lead to the appearance of causal pathologies, i.e. all 
future-directed world lines remain unclosed. In the sec- 
ond > 1, the situation is cardinally changed (this 
case is illustrated by the spacetime diagram in the figure 
1). 




Fig.l. 



Namely, now there is such a moment of time (the value 
of 77 = 77*) when one of the cone's side takes a 'horizontal' 
position, i.e. one of the null curve's equations becomes 
77 = const. But these lines are closed (see the Eq.(Q)), and 
hence, at the moment 77* the closed null curves appear in 
our model. We shall speak that a time machine is being 
formed at this moment of time. If 77* — 00 (a = 1) 
the time machine is formed in the infinitely far future. 
Otherwise, it is formed at the time 77* < 00. In this case 
at later times rj > 77* the closed line 77 — const lies inside 
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of the light cone, i.e. the region 77 > 77^ contains closed 
timelike curves. 

Thus we may conclude we have constructed the model 
of spacetime in which, depending on the parameter a , 
the chronology horizon is being formed at some moment 
of time. 

Now let us investigate a behavior of quantized fields in 
this model. 



III. A PARTICLE CREATION 

A. Scalar field: solutions of the wave equation 

Consider a conformal massless scalar field for which 
the Lagrangian is 



(11) 



The scalar field obeys the wave equation 

□ = 0. (12) 

In addition, it follows from the identification rule (|t]) and 
from the quadric form of the lagrangian ( |TT| ) that the 
scalar field has to obey the periodic condition (the ordi- 
nary field) 

0(r7,£ + L) = 0(r7,£) (13) 
or the antiperiodic one (the twisted field) 

+ L) = -<j>( V ,0- (14) 
In the metric (Q) the wave equation ( |l2| ) reads 

(1 - a 2 )dl + 2ad v d i - S| - 2aa'd v + a'd^l 4>(t], f) = 0, 

(15) 

where 8^ = ^,8^ = -^ and a prime denotes the deriva- 
tive on 77, a' = da/drj. 

First of all, let us solve this equation in the asymp- 
totical regions. In the in-region, where a(rj) — > and 
a'(r]) — > 0, the Eq.(^5|) reduces to 



- df Mr?, £) = 0, (16) 

The complete set of solutions of this equation is 

0(±,in) = £(m) e ;fc„c eT ^ 

where 



27T77, 

k n = —jr~, n = ±1, ±2, 



(ordinary field) 



2tt(71 + i) 

— - - 2 \ n = 0,±1,±2,. 



lu = \k n \ and 0i + ' m ^ and 0i~' m ' ) are positive and negative 
frequency solutions (the "in- modes" ) in the in-region, re- 
spectively. The in-modes form the basis in Hilbert space 
TL with the scalar product 



T>n , <Pn' , 



indued!?, 



(18) 



where c?E M = dS n^, with dY, being the volume element in 
a given spacelike hypersurface, and being the timelike 
unit vector normal to this hypersurface. Choosing the 
hypersurface 77 — const we may write down the scalar 
product (|T|) in the in-region as follows: 



') = -*/ d£ 
Jo 



9?7 dij 1 
The in-modes are orthonormal provided 

£)(in) = 



77— const 



1 



■\/47r|n + a\ 



(19) 



(20) 



where a = for an ordinary scalar field and a = \ for a 
twisted one. 

An alog ously we may find a solution of the wave equa- 
tion ( Jig ) in the out-region. There a(n) — > a and 
a'(?7) — » 0, and the Eq.(|i~5|) reduces to 



(1 - a 2 )d* + 2a dr,dt - 9|j 0(77, £) = 0, 
The complete set of solutions of this equation is 

^(±,out) _ £)(out) e ifc„5 e -ifc„a /3 _1 i) e ^iw/3 _1 r) 

where we denote 

= 1- al 



(21) 



(22) 



(23) 



and 0l +,out ^ and (f>n ,out ^ are positive and negative fre- 
quency solutions (the "out-modes") in the out-region, re- 
spectively. The out-modes ( |22] ) form the basis in Hilbert 
space Ti. with the scalar product ( |l8|) which in the out- 
region reads 











dr, Or, ™ 



l-di 



®£ ®£ /J r/=const 

The out-modes are orthonormal provided 



£)(out) 



Ait\n + a\ 



(24) 



(25) 



(twisted field), 



Now let us solve the Eq.([L5|) in a general case. Noting 
that the metric (Q) is invariant under translations in the 
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^-direction and takin g in to account the periodic or an- 
tiperiodic conditions (|l3|), ( [hi] ) we may find solutions in 
the following form: 



<M?7,£) = u n (r))e 



(26) 



Substituting the expression (|2q ) into the wave equation 
( |l5| ) we obtain the equation for u n (rf): 

(1 - a 2 )u" + 2a(ik n - a')u' - ik n {ik n - a')u = 0. (27) 

Introduce a new function v(r)) by the relation 



u(rj) = v(r/) exp 



71 a(ik n - a') 
l-o 2 d " 



(28) 



Note that in the in-region the relation (J2§|) has the 
asymptotical form 



u(rj) = v(r)), 
whereas in the out-region it reads 

u(v) = v{r 1 )e- lk ^ a °^ 1, i 



(29) 



(30) 



After substituting Eq.(28|) into Eq.(|27|) we obtain the 
second-order differential equation for 11(77) in a so-called 
normal form v" + Vl 2 (rf)v = 0: 



k 2 n + a 12 + a"a(l - a 2 ) 



(1-a 2 ) 



2\2 



v = 0. 



(31) 



Further we shall solve this equation only for modes for 
which the conditions 



K » 



k 2 n » a"a(l - a 2 ) 



(32) 



are fulfilled. That is we shall only consider the modes 
whose wave length is much less than a typical scale of 
variation of the function 0(77). Taking into account the 
conditions (^) we could now neglect in ( |3l"| ) the terms 
a' 2 and a" a(\ — a 2 ) and rewrite^] 



k 2 v 



(I -a 2 ) 



2\2 



0. 



(33) 



Now let us restrict our consideration by the special form 
of the function a(rj): 



a2 (v) = ^o(l + tanh 7 77), 



(34) 



where 7 is a parameter. It is easy to see that the function 
a(rj) defined by Eq.(|34|) possesses the necessary asymp- 
totical behavior (0). The quantity T ~ (27) ~ 1 gives the 



*Let me emphasize that the approximation which we use 
is no WKB approximation. Remind that the last means, 
roughly speaking(!), neglecting terms with Q' , Q" , ... in the 
solution of the equation v" + Q?{rr)v = 0. 



typical time of variation of the function 0,(77) from one 
asymptotical value to another one. The conditions (|3^ ) 
reduce now to 



» 



27 T 2 ' 



(35) 



Substituting the expression ( |34| ) into Eq. (|33j) we could 
rewrite the equation ([33|) as 



4fc 2 w 



(2 - a 2 (l + tanh7ry)) 2 



0. 



(36) 



A solution of this equation could be found in terms of 
hypergeometric series. In paticular, to find a solution 
which will be positive frequency in the in-region we have 
to choose the solution of the equation (B6|) as follows: 



V^(V) = DtH-z)-^(l-pzrF(q,r;s;pz), (37) 



where d!" 1 ' is defined by 




(38) 



and F(q,r; s; z) is a Gaussian hypergeometric function. 
Taking into account the relations (|2^) and ( p6| ) we can 
now write down the solution of the wave equation jl5| ) 

as 



x(l-/3zyF(q,r;s;l3z) 



71 a{ik„ 



-dfj 



(39) 



Note that F(q, r; s; 0) = 1 for any q, r and s. Now it is not 
difficult to see that this solution in the in-region, where 
77 — > — co or z — > —0, has the following asymptotical 
form: f n (iliC) ~ Dn * e lkn ^ e - lujr i ( which coincides with 



the positive frequency in-modes 



B. Bogolubov coefficients 

Now let us remind some mathematical aspects for de- 
scribing of the physical phenomenon of partical creation 
by a time-depend gravitational field. 

So, we have obtained the set of solutions {f n } which 
are positive frequency (the "in-modes") in the past. Let 
{F n } be positive frequency solutions (the "out-modes") 
in the future. (We do not need to know an explicit form 
of these solutions. It will be enough to know their asymp- 
totic properties in the out-region, i.e. F n w 0n +,out ^.) We 
may choose these two sets of solutions to be orthonormal, 
so that 
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(f n J r :,) = (F n ,F*,)=Q. 



(40) 



The in-modes may be expanded in terms of the out- 
modes: 



(41) 



Inserting this expansion into the orthogonality relations, 
Eq. (pip]), leads to the conditions 



and 



(42) 



(43) 



The field operator, <f>, may be expanded in terms of either 
the {/„} or the {F n }: 

<P = VW" + 4 /*) = Y\{KF n + blF*). (44) 



The a n and a) n are annihilation and creation operators, 
respectively, in the in-region, whereas the b n and b^ n are 
the corresponding operators for the out-region. The in- 
vacuum state is defined by a n |0)j n = 0, Vn, and describes 
the situation when no particles are present initially. The 
out- vacuum state is defined by 6„|0) OU f = 0, Vn, and 
describes the situation when no particles are present at 
late times. Noting that a n = (</>, /„) and b n = {4>,F n ), 
we may expand the two sets of creation and annihilation 
operator in terms of one another as 



■8* b ] ) 



or 



E 



(45) 



(46) 



This is a Bogolubov transformation, and the a nm and 
B nm are called the Bogolubov coefficients. 

Let us assume that no particle were present before the 
gravitational field is turned on. In the Heisenberg ap- 
proach |0} !n is the state of the system for all time. How- 
ever, the physical number operator which counts parti- 
cles in the out-region is N m — b^bm. Thus the number 
of particles created into mode m is 



(N m ) = in (0|6j n 6 m |0) i 



(47) 



To find the Bogolubov coefficients in our case we have 
to determine an asymptotical form of the in-modes in the 
out-region where i) ^ oo or z -> - oo. With this aim we 



consider the analitical continuation of the hypergeomet- 
ric function F(q, r; s; z) into the region of large values of 

1*1 @ 

F(q,r;s;z) = 

1 (r)l (s — q) z 

+ W0) { - zrrF ^ i - s+r]i - q+r] - z y (48) 

After substituting this expression into the Eq.([59j) it is 
not difficult to see that the asymptotic of the in-modes 
in the out-region is 

/«(»?, « D<g^j k ^e- ik " a °P~^ 



(-/?)" 



pi/2T(r)r(s - q) 



-iuj/3 77 



-0) 



r( S )T(q-r) 
P l / a T(q)T(8 - r)' 



Ju0 Tj 



(49) 



Comparing the last expression with the asymptotical 
form of the out-modes, Eq.(J23), we may write 



= A n F n + B n F*_ n 



where 



A r , 



i-0) 



/ 3 1 /2r(r)r( s - 9 ) : 



Bn = ( py U 8-i(i+8) r(a)r(g-r) 



/3 1 /2r( g )r( s - r ) ' 



(50) 
(51) 
(52) 



As follows from Eq.([4l|), the coefficients A n and B n 
are the related to the Bogolubov coefficients by a nm = 
A„S nm , Pnm = B n 5 n ,- m - The number of particles cre- 
ated into mode n is now determined as 



(N n ) = \B n \ 2 . 
Using the Eq.(|52|) gives 

T( S )T(q 



(N n ) = 3- 



T(q)T(s - r) 



(53) 



(54) 



To carry out calculations in the above expression we shall 
use the following formulae |13] : 



y sinh ny ' 



\rC 1 + i y)\ 2 = 



cosh Try ' 



|r(i + iy)\ 2 = - 



sinh Tty 



(55) 



By using this formulae and the relations ( J3q ) we may 
finally obtain the spectrum of created particles: 



(Nn) = 



7/3 



COSh7T-l/( 



7/3 , 



- 1 



2 sinh ^ sinh ^ 

7 7/3 



(56) 
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C. Particle creation near the chronology horizon 
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Now let us analyse how much particles of the massless 
scalar field is created at the moment of the chronology 
horizon forming. As was mentioned above, the value of 
the parameter a determines either the time machine is 
formed or not. If a < 1 then the chronology horizon is 
absent in the spacetime, but if a a = 1 then it appears 
in the infinitely far future. The expression (||j) gives us 
the spectrum of particles created in the out-region. The 
spectrum depends on the parameter a (/3 = 1 — and 
in the limit a — > 1 the expresion ([56]) will determine the 
number of particles created near the chronology horizon. 
Going to the limit a a — > 1 in Eq.d5^) gives 



1 



sinh ' 

7 



(57) 



Thus we see that the number of particles created into 
mode n near the chronology horizon is finite. We may 
also conclude that the full number of particles N — 
J2 n (N n ) will be finite because the spectrum ( ^7| ) is ex- 
ponentially decreasing. 



IV. CONCLUSION 

Let us summarize. In this work we have constructed 
the model of spacetime in which, depending on the 
model's parameter a , the chronology horizon could be 
formed; there are no causal pathologies if a < 1, and 
the chronology horizon appears at the moment of time 
?7* < oo if a > 1. In the case a Q — 1 closed lightlikc 
curves are formed in the infinitely far future. The model 
represents a two-dimensional curved spacetime with the 
topology R 1 x S 1 which is asymptotically flat in the past 
and in the future, but which is non-flat in the intermedi- 
ate region. As a consequence, in this spacetime the cre- 
ation of particles of quantized fields by the gravitational 
field is possible. We have studied the particle creation 
of a massless scalar field in the case a a < 1, i.e. in the 
case when the spacetime is globally hyperbolic and has 
no causal pathologies. As a result, the spectrum of cre- 
ated particles has been obtained (see Eq. ([Scil) ) . In the 
limit a — > 1 this spectrum gives the number of par- 
ticles created into mode n near the chronology horizon 
(see Eq.(p7|)). The main result we have obtained is that 
the number of scalar particles created into each mode 
as well as the full number of particles remain finite at 
the moment of forming of the chronology horizon. This 
result might mean that the phenomenon of particle cre- 
ation could not prevent the formation of time machine. 
However, to do the final conclusion one has to take into 
account a backreaction of created particles on a space- 
time metric. 
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